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Abstract. A general quantum adiabatic theorem with and without the time- 
dependent orthogonalization is proven, which can be apphed to understand 
the origin of activation energies in chemical reactions. Further proofs are also 
developed for the oscillating Schwinger Hamiltonian to establish the relation- 
ship between the internal (due to time-dependent eigenfunctions) and external 
(due to time-dependent Hamiltonian) timescales. We prove that this relation- 
ship needs to be taken as an independent condition for the quantum adiabatic 
approximation. We give four examples, including logical expositions based on 
the spin- 1 two-level system to address the gapped and gapless (due to energy 
level crossings) systems, as well as to understand how does this theorem allows 
one to study dynamical systems such as chemical reactions. 
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1. Introduction 

The quantum adiabatic theorem (QAT) and its approximation (QAA) have been 
the backbone in many areas of quantum physics, namely, in condensed matter 
theory (via the Born-Oppenheimer approximation) in atoms, molecules and 
quantum chemistry in quantum field theory via the Gell-Mann and Low the- 
orem [3], and presently in the adiabatic quantum computation The theorem 
was first discussed by Ehrenfest [5 , and later was formally derived by Born and 
Fock [B] and Kato [Tj. Other modern proofs for QAT can be found in Refs. [51111 
Uni HH Ha miUj. We refer the readers to Refs. HH [H] for thorough reviews 
on the development of QAT. Apart from the historical perspectives, Comparat [11] 
have also proven the existence of two sufficient QAT conditions for the oscillat- 
ing Schwinger Hamiltonian. Teufel-Spohn [TJl [TB] and Avron-Elgart [T^j on the 
other hand, have provided the proof for the validity of QAT for systems obeying 
the Born-Oppenheimer approximation, and for gapless systems, respectively. We 
will revisit the two conditions given by Comparat to prove that these conditions 
are in fact equivalent to the relationship between Tex and Tin. Here, Tex and Ti„ 
are the respective external and internal timescales. Subsequently, we will discuss 
how the new theorem proven here addresses the gapless system. Due to numerous 
applications of QAT in physics, a proper theorem with explicitly defined conditions 
is needed. Here, we prove the existence of such a theorem, which can also be ex- 
ploited to understand the origin of chemical reactions. To arrive at the standard 
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QAT for gapped system, we first let H{t), tlj{t) and E{t) to be the time-dependent 
Hamiltonian, eigenfunction and eigenvalue, respectively, such that [5] 

(1.1) H{t)^^it) = Ejt)Mt), 

and {ipnlt)} constitutes a complete orthonormal set that satisfies 

in which ^n{t) can be expressed as a linear combination 



(1.3) 

with the phase factor 
(1.4) 

and {n, m} G N*. We use N* and M to denote the set of positive integers excluding 



9„(<) = -^^ E^{t')dt\ 



zero and the set of real numbers, respectively. Here, Eq. ( 1.3 ) is the general solution 
to the time-dependent Schrodinger equation 



(1.5) 



Therefore, one can show that if the initial Hamiltonian, i?init(i) evolves very slowly 
to a new final Hamiltonian, Hfin{t), then a particle that starts out in the n^^ 
eigenstate (*L""(0 = V'«(i)e'^"(*)) of iJi„it(t) is expected to occupy the ground 
state (^'^"(0 = '0„(i)e*^"(*)e*^"(*)) of i7fi„(t). Observe that the final eigenstate has 
non-trivially, picked up only a harmless phase factor 



(1.6) 



In the above formulation, Griffith used Eqs. (1.3) and (1.5) to obtain the required 
coefficient (in explicit form) j9j. 



Cm{t) = -Crn{t){iprn{t)\iPm{t)) - ^ C„(t) 



iprn{t) H{t) (p.a{t) 



,»(e„(t)-e™(t)) 



En{t) — Ejnit) 



(1.7) 



that can be used to find out how slow the evolution should be, so that the initial 
particle still occupies the ground state of the new (final) Hamiltonian. Here, Cm(t) 
and Cn{t) denote the m*^ and n*^ eigenstate coefficients, respectively, such that 



|cm(t)P + |cn(OP = 1- From Eq. (1.7), we can notice that both En{t) and Ejn{t) 
exist at all times, and their eigenf unctions are orthogonal. This means that if En{t) 
and Em{t) are degenerate at certain point of time, then Eq. (1.7) cannot be used 
to derive the adiabatic criterion. The reason is that the gap that exists when t = to 
does not exist for t — ti. Now, assuming non-degeneracy, the QAT criterion is 
given by 



(1.8) 



(p„i{t) H{t) (pnit) 
En{t) — Em{t) 



< Cmit), 
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which needs to be satisfied for one to exploit QAA. Here, fm{t) is the usual t- 
dependent eigenfunction for the m'^ eigenstate and £'„ {t) is the t-dependent eigen- 
value for the n'^ eigenstate, and the dot represents time derivative. There are two 
important assumptions associated to the original QAT given in Eq. (1.8 1. First, 
there is no quantum phase transitions or chemical reactions when the initial quan- 
tum system evolves from to ti. Secondly, the QAT is only applicable for a two- 
or a multi-level system with well-defined energy gaps (strictly no energy level cross- 
ings). For example, when one considers a two-level system for simplicity, then there 
exist two i-dependent eigenf unctions, observable at all times, and they are always 
nondegenerate and the two eigenvalues are En(t) ^ Em{t). In fact, we can add as 
many eigenf unctions as we possibly could (as denoted by n and m), but we chose 
two, to keep the discussion straightforward. Of course, later we will go through the 
examples with large number of energy levels and energy level crossings. 

For the sake of argument however, if we assume that En{t — to) ^ En{t — ti) 
(this is definitely different from En{t = t) ^ Em{t = t), as described earlier) then 
we cannot apply the quantum adiabatic approximation [given in Eq. (1.8)] simply 
because both of these eigenvalues {En{t — to) and £'„(t = ii)) are not observable 
simultaneously. Recall that to < ti. Furthermore, En{t = tg) ^ E„(t = ti) implies 
the energy gap (g) in the form of .g = \En{t = to) — En{t = ti)\. The energy gap in 
this form cannot be defined as an energy gap because 

En{t = to) — En{t = ti) 

h 

E.,{t)^E^it) 



(1.9) 



Jit) = 

'^(i)Bolir 



1 



In other words, physically valid energy gaps must satisfy w(t)Bohr- Later, we will 
show that uj'{t) can be associated to the possibility of a particle traveling forward 
and backward in time. Now, it is also important to note here that if the external 
disturbances with a timescale, 



(1.10) 



hr 



Tox(t) 



CX (ipmit) H{t) (Pn{t) 



approaches WBohr due to T^^ ~ Tim then one cannot apply QAA. Now, one should 
not be carried away into thinking that Eq. (1.8) lacks a small parameter, which is 
required to understand the quantum adiabatic evolution. For example, the small- 
ness is usually controlled by the dilation operator, r such that Sex = Tq^/t. In our 
formalism however, this smallness is controlled by the relationship between Tin and 
Tox, without the need to introduce any ad-hoc control parameter hy hand. We will 
get to this point in detail in Example 1 when we proof Eq. (1.10). The adiabatic 
criteria, Tex 3> Tin and Eq. (1.8) have been proven in Ref. [5], however, without the 
proper definitions for these timescales. The existence of uj' {t) means that we are 
equating ipn(t) with ipn(t = to) and (pmit) with ipn{t = ti) in Eq. (1.8). In this case, 
we are referring to only one eigenfunction at all times such that (^„(i = to) evolves 
to (fnit = ti), where (pnit = to) and (pnit = ti) may or may not be orthogonal to 
each other. If they are orthogonal, then one obtains u!'{t). On the other hand, if 
they are not orthogonal, then uj'{t) does not exist. 

In addition to our arguments based on Eq. (1.9), Comparat |llj have also shown 
that the condition given in Eq. (1.8) is valid provided that the Hamiltonian is real, 
gapped and non-oscillating. However, the generalized conditions derived there [TT] . 
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which are vahd for any Hamiltonian, do not consider gapless systems nor the Hamil- 
tonians that allow t-dependent orthogonalization. In addition, the internal and 
external timescales are also not properly taken into account. The i-dependent 
orthogonalization and large changes to the Hamiltonians (due to large external 
disturbances) are essential because they can lead to continuous or second-order 
quantum phase transitions (QPT) during chemical reactions. The readers are re- 
ferred to Refs. [101 mi m 113] for further details on the concept of QPT and its 
applications in strongly correlated matter. 

We anticipate that this continuous QPT proposed by Sachdev [IHl [H] and 
Senthil [HI US] can also occur due to the existence of energy level crossings due 
to large external disturbances, with or without the t-dependent orthogonalization, 
which could be the origin of any chemical reaction between two chemical species. 
This means that the energy levels of unreacted species must cross by overcoming the 
activation energies to form new compounds, and these new compounds (if formed) 
will have a new set of energy levels or bands. If the formed compounds are of 
non free-electronic systems, then there is a theorem based on the ionization energy 
theory that states— the energy level spacing of such a system is proportional to 
their constituent atomic energy level spacing [24^. Physically, this means that the 
energy level spacing of this newly formed compound is proportional to the energy 
level spacing of their constituent chemical elements. 

Before we move on, let us first try to understand the concept of internal and 
external timescales in the QAT. It is well known now that there are different versions 
of QATs available for different applications as reviewed by Teufel [5], and Avron- 
Elgart [H]. However, all the existing QATs are solely based on the scaled time, 
Sox = ^cx/tox, Sox € [0, 1], without explicitly referring to the internal and external 
timescales. In particular, Scx here is independent of the internal timescale. The 
reason for this scenario is that there are no available equations to properly define 
Tox and Tin independently. Anyway, the scaled (dimensionless) time, Sox contains 
the dilation time or a time dilation operator, Tcx, such that r^x allows us to slow 
down the time tex to obtain a slowly varying Hamiltonian with respect to the value 
1 (when Tex = iex)- Hence, tox is the original (unsealed) physical time for the 
Hamiltonian, and t^x determines the dilation of this external physical time. 

In the unsealed formalism, we simply write Tox to denote external timescale, to 
avoid confusion with t-dependence in H{t) and ifi(t) (eigenfunction) . Otherwise, we 
will have to write H(Tcx) and (p(Tin), and this notation will confuse the readers into 
thinking that both H{Tcx) and 95 (Tin) are not normalized to the real-space physical 
time, t, or V'(7in) is from a different Hamiltonian, and so on. As a consequence, we 
can see that the Hamiltonian with scaled time, H{scx) does not take the internal 
timescale into account, which can only be taken into account by defining the time 
variation in the wave functions or the eigenfunctions for a given H{scx)- Similar to 
Sox, we can also define Si„ — tin/Tin, where Sm S [0,1]. In these scaled notations, 
tin refers to the original physical time of the wave functions, and Tin determines 
the dilation of the internal physical time. Here, Sin is independent of the exter- 
nal timescale. Again, in the unsealed formalism, we write Tin to denote internal 
timescale. 

Having understood the different notations used in the scaled and unsealed for- 
malisms, we can now properly define the time, Tox as the external timescale (char- 
acteristic time for changes in the Hamiltonian), while Tin is the internal timescale 
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(characteristic time for changes in the wave function or eigenfunction) . In order to 
understand the existence of these properly defined timescales, we do the following: 
we can define = tcx = t and write Sin = i/rin = t/Tln and Sox = t/rcx = t/Tcx- 
This normalization transforms H{so-x) back to H{t), and the only thing left to do 
here is to find the inequality relation between Tex and Tin. However, this _ff (scx) or 
H{t) still ignores the internal timescale (because it is fixed as a constant), which 
leads to an important question— how small should (or T^^) be with respect to 
Tin in order to violate the QAA? This question means that we cannot claim a given 
Hamiltonian varies slowly if — >■ oo because we have to make sure Tin does not 
go to oo when Tcx — >■ oo. We will prove here that the relation between Tin and Tex 
needs to be treated as a separate condition altogether. 

Now we will discuss why the concept of quantum adiabaticity is correct, and why 
any violation that may exist only indicate the non-applicability of the QAA for cer- 
tain systems. The existence of non-orthogonality in strongly interacting physical 
systems (or non-linear quantum systems) has been discussed by Yukalov |17| as 
a result of non-linear terms in a given Hamiltonian that destroys the adiabatic 
evolution. These non-linear terms can originate from various interactions such as 
electron-electron, electron-phonon, spin-orbit, spin-spin and so on, which may give 
rise to numerous energy-level crossings along the time line, t. As anticipated, this 
means that we need to keep track of the orthogonality of the eigenfunctions with 
respect to time (which will be explained in the subsequent sub-section) . It is to be 
noted here that strongly correlated systems, including quantum dots and nanos- 
tructures (non-free-electron systems) do not form gapless systems. For example, 
the electrons and the phonons cannot be decoupled, unless Tdcctron ^ Tphonon [T5] . 

Any violation to a given QAT means the existence of a non-adiabatic process, 
and this does not imply the QAT is incorrect. Violating the QAT deliberately is 
important as a crosscheck because if the QAT is violated for adiabatic processes, 
then the QAT is either (a) not applicable for that process (for example, applying 
gapless or gapped QAT to a system that has time(i)-dependent orthogonalization), 
or (b) the said QAT is insufficient (again, this does not mean the concept of quantum 
adiabaticity is incorrect). Consequently, deliberately violating a given QAT can be 
used as a mean to crosscheck the correctness and sufficiency of the conditions given 
in a particular QAT, as well as to have a deeper understanding of the conditions 
needed to guarantee the QAA, for example, with respect to the existence of self- 
orthogonality [TO] . 

Here, we develop a new version of the QAT that (i) allows t-dependent orthog- 
onalization for each single eigenfunction as an additional condition, which can be 
used to explain the possibility of a particle traveling backward in time, (ii) allows to 
understand the origin of chemical reactions in the absence of, and in the presence 
of i-dependent orthogonalization, and (iii) gives a proper definition for the Tex and 
Tin and establishes the relationship between them as a valid and separate condi- 
tion. Points (i) to (iii) are the main results of this paper. In addition, we apply 
our theorem to degenerate energy levels due to energy level crossings to validate 
the QAA with properly defined Tex and Tin without invoking any Rabi oscillation. 
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2. Generalized Quantum Adiabatic Theorem 



We use "H for the complex Hilbert space, while \ \(p\\ is the norm of an eigenfmic- 
tion {(p). 

Theorem 1. The quantum eigenf unctions (or wave functions) that satisfy the 
quantum mechanical postulates 25 are represented by the orthonormalized com- 
plex vectors, \ip{t)) = (•••, W(tj)),- ■ ■ , \(p{tk)), ■■■ , l^fiitz))), where \ip{t)) £ H and 
||(y5(t)|P = {Lp{t)\tf{t)) = Jjj^ =1. If a quantum system with a time- 

dependent Hamiltonian, -ffinit(^) is initially {t ~ tj) in its ground state {\ipn{t)) or 
\^m(tj))), evolves into a final {t — tk) quantum system with a Hamiltonian, -fffin(t) 
then the probability for the final quantum system to be in its ground state {\ipm{t)) 
or \ipmitk))) is controlled by these three criteria 



(2.1) 



Vra{tk) H ipmitj)^ - Eyn{tj){ip„i{tk)\(prn{tj)) 



Em{tj) — Ejnitk) 



(2.2) 



¥>„(*) H{t) ifnit) 



En{t) — E„i{t) 



< Crn{t), and 



Tcx(t) > Tin(t). 



(2.3) 

Here t S [tj,tk] G K, and {^Pm{tk)\y^m{tj)) = ^t^tj ~ implies ipm{tj) is orthogonal 
to ipm{tk)- On the other hand, {(pm{tk)\^m{tj)) — ^t^tj — 1 implies ipm(tj) is not 
orthogonal to ipm{tk) and {j, fc,r7i,n} € N*. 

Remark 1: Equation (2.1) keeps track of all the individual t-dependent eigen- 
functions in a given system. QAA is strictly valid if (a): Cm(^fc) in Eq. (2.1) is 
replaced with = Cm{tk), (b): < Cm{t) in Eq. (2.2) is replaced with <C Cm{t), and (c): 



> in Eq. (2.3) is replaced with In fact, we can reversibly switch to strict QAA 



conditions if needed, or loosen them, or use them to invalidate the QAA in certain 
systems. We can now write a statement in regards to QAT: The probability for a 
particle that was initially in the ground state of an initial Hamiltonian, -ffinit(Oj to 
occupy any excited state of some final Hamiltonian, Hfi^{t) is approximately zero, if 
^?init(i) changes slowly to its final form, 7Jfin(i) such that Eqs. (2.1), (2.2| and (2.3) 
are satisfied. 

Remark 2: Time discretizations have been carried out to obtain instantaneous 
eigenfunctions. In other words, the t-dependent eigenfunction has been time- 
discretized in the form of fm{t) — ■ ■ ■ , fmiij), ■ • • , ^m{tk) ■ • • , ^m{tz)- Time dis- 
cretizations are carried out in pairs {tj and t^) to obtain two time-discretized eigen- 
functions (ipm{tj) and ifmitk)) so as to check their orthogonality. The sole purpose 
of this time-discretization is to capture the t-dependent orthogonalization and/or 
the wave function transformation (that will be exposed in the last section). 

Remark 3: In addition, it is to be noted here that c{tj) = c{t = tj), f{tj) — (p{t 



tj) and 9{tj) — 9{t — tj). As for the time derivative variables, c{tj 
ip{tj) = ^Q^*-* (t = tj), and so on. 



_ 9c(t) / , _ 



dt 



it = h), 



Proof for Eq. (2.1): The validity of Eqs. (2.2) and (2.3) have been proven in 
Ref. The discussion on how Eq. (2.1) is related to chemical reactions is given 
in Example 4, while the proof for Eq. (1.10 1 is given in Example 1. We will start 
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with a single eigenfunction (to = 1) for two discretized times, i.e., when t — tj and 
t = t^. For simpUcity, we let j — Q and fc = 1 from here onwards, and then we 
can just repeat this procedure for j = 1 and fc = 2, j = 2 and fc = 3, and so on. 
Moreover, (pm{to) is allowed to evolve and to be orthogonal or not orthogonal to 
Vm{ti) when t — tl. The ^-dependent Hamiltonian, say for t — to is given by 



(2.4) 



HitoMto) - EitoMto). 



The i-dependent Schrodinger equation is given in Eq. (1.3) in which, ^'(i) can be 
expressed as a linear combination of ipito) and (p{ti) (from Eq. (1.5)) 



(2.5) 



^{t)^J2<*oM^o)e'"''°^ ^ te[to,ti]. 



The summation is used between smd ti because we have discretized the t- 
dependent eigenfunction (See Remark 2). Here, we also have c(io), which is a 
coefficient that is related to the probability of finding a particle at a given point of 
time (io or ti), and with a particular eigenvalue, Em{to) or Em{ti), respectively, for 
TO = lor 2 or • • • . For example, Eq. (2.1) implies |cm(io)P + |cm(ii)P — 1, while 
Eq. (2.2) implies |cm(i)p + |c„(i)p = 1 (for to ^ n). Recall here that we are only 
considering a single eigenfunction (to = 1), and due to Remark 2, it is convenient 
for us to drop the subscript to. From Eqs. (1.3) and (2.51 

tl 

* = E HtoMto) + c{to)ip{to) + ic{toMto)e{to)] e^'^'^K 

to 

tl 

ihJ2 [KtoMto) + c{to)ipito) + ic{toMto)0ito)] e*''(*°' 



(2.6) 

Using 

(2.7) 

we obtain 



J2(^ito)Hip{to) 



Oito) 



E{t')dt' 



0{to) 



:E{to), 



iHj2 [c(to)(^(io) + c(io)¥'(io)]e*'^*°^ 



to 



J2cito)EitoMto)e^"''^'^ 



(2.8) 



to ^ 

tl 

ihJ2 [citoMto) + citoMto)] e''^'"^ 

to 

tl tl 

+ J2 c{to)E{toMto)e^'^^'°^ = J2 cito)EitoWito)e'"^'"^ 



to 
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which leads to 

(2.9) |^c(io)¥'(io)e^'(*"^ 



ie(to) 



^0 *0 

Now comes the crucial part, we will take the inner product with (p(ti), which is an 
eigenfunction at a later time, ti (ti > to), evolved from <p(io)- Therefore, 

(2.10) {v{ti)\J2 c(io)|^(io))e''(*«) = -{^{h)\ f2 c{tomto))e''^'''\ 



to 



where 
(2.11) 



From Eq. (2.11 1, one obtains 

(2.12) 1^ i{to)6t,toe''(">^ =-J2 <to) ivitimh)) e^'^'-\ 

to to 

Now, let us assume that the evolved eigenfunction, ip(ti) is not orthogonal to if(to) 
or St-^to = 1- Hence, Eq. (2.12) can be written as 

c(ii)5t,t„e'«(*^) = - Ec(to) {^{hmto)) e'^'(*°) 

to 

tl 

(2.13) c{h) = -J2cito){^ih)mo))e^^'^'»^-'^'^^l 



Differentiating Eq. (2.4 1, taking the inner product withi^(ii) and using {ip(ti)\H\ip(tQ)) 
E{ti){ip{ti)\(p(tQ)) we can derive 

H^{to) + H^{to) = E{toMto) + Eito)^{to), 

(2.14) {ipiti)\H\^ito)) + Eih){ifih)\ip{to)) = Eito) + E{t^){^{h)\^{to)). 

Since the single eigenfunction at to remains as a single eigenfunction at ii, there 
is no available excited eigenfunctions for any transition to occur. Therefore the 
transition probability is simply zero. Invoking 5t-^^to = 1j can rewrite Eq. (2.14) 
as 

(2.15) {E{to) - E{t,)){ip{ti)\ip{tQ)) = {ipiti)\H\cpito)) ~ E{to). 
Substituting Eq. (2.15) into Eq. ( 2.13[ ), we obtain 

c{t,) = -c{t,){^{t,)m,)}- ^ c(io)(^(ii)|^(io))e'[''(*°)-''(*^)l 



(2.16) =-c(ti)(^(ti)|^(ti))- ^ c{to) 

to^ti 



H 



^{to))-E{to) 



,i[e{to)-B(ti)] 



E{to)-E{ti) 



Equation (2.16) is straightforward where c(to) = c(ii) — 1 because there is only 
one eigenfunction, and is not orthogonal along the time line for t e [io,^i]- Let us 
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now invoke the orthogonality {Stj^tg = 0)j which imphes that (from Eq. (2.13)) 

tl 



(2.17) 



,ie(to) 



Using Eqs. ( |2.17| ) and ( |2.14[ ) 
(2.18) cih) = - 



(^(ii)|^(<i)) ^ ii;(to)-s(ii) ' 

which in turn allows us to conclude that 



(2.19) 



^ih) H ifito) 
E{to) - E{ti) 



1 ^ c{ti) « 1, or 



(2.20) 





H 


f{to)) 




- E{h) 



< 1 =^ c{ti) < 1. 



Equations ( |2.19 1 and (2.20) are equal to the adiabatic criterion given in Eq. (2.1) 
ifi?(to)(^(ii)|^(io)) v^OB 

As indicated earlier, systems with a single eigenfunction, and in the presence 
of ^-dependent orthogonalization, we will face the consequence of a particle trav- 
eling backward in time due to c{ti) <C 1 (Eq. ( |2.20[ )) in accordance with w'(t) 
(see Eq. (1.9)). This is physically not possible because we will not be able to or- 
thogonalize (assuming this can be done) the eigenfunction without disturbing the 
occupied electron. In other words, it is not possible because the eigenfunction here 
itself represents the wave function of the electron based on the Hermitian quantum 
mechanics [23]. Therefore, Eq. (2.20) strictly implies that the ground state electron 
has been excited or ionized. We cannot be sure how this (excitation or ionization) 



is possible from the criterion given in Eq. (2.1 ) alone, it could be due to degeneracy 



or invalid Eq. (2.3) [for example T^xit) < T'in(i)]- Later in Example 4 we will justify 
that this phenomenon (c(ii) <^ 1) can be related to chemical reactions (with and 
without the ^-dependent orthogonalization). 



3. Applications of Theorem 1 

Now we are ready to discuss the applications of Theorem 1 , and how this theorem 
relates to these examples both mathematically and physically. The first example 
is for spin-1/2 two-level system, which is by far the most straightforward one. The 
second one deals with the same spin-1/2 two-level system but with more and compli- 
cated <-dependent external disturbances. Subsequently, we discuss gapless system 
due to energy level crossings, while the last one invokes the t-dependent orthogonal- 



ization in the spin-1/2 two-level system via Eq. (2.1) to check the applicability of 
QAA for such systems. We also show that Eq. ( 2.1 ) can be exploited to understand 



chemical reactions in the absence of t-dependent orthogonalization. 
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Figure 1. A constant magnetic field, B (magnitude Bq) rotates at 
a constant velocity, il, and with a constant opening angle, a. Here, 
B(t) = BQ[smacos{ilt)i + sin a sin(rit)j + cosak]. See Example 1 
for details. 



Example 1: Spin-i two-level system in the presence of magnetic field. 

The first example is the one derived by Griffith [9 . An electron can have spin 
up or down, and therefore, in the presence of rotating (constant velocity = f2) 
and constant magnetic field with a constant opening angle, a, one can control the 
transition probability of this electron (from spin up to down) to satisfy |9j 



(3.1) 



|(*(i)l^dow„(i))l' 



where A^ = il^ + io^ohi- ^ Sri^Bohr cos a and i:JBohr = [Enp - i?down]/^ = l/Tln- 
Obviously, in this example we have f2 = l/T^x- In the limit Tox ^ due to f2 <C 
'^Bohr: one obtains the transition probability to spin down to be zero because A — >■ 
'^Bohr- Therefore, (^i/A) — )• (O/tLiBohr) — in accordance with Eq. (3.1), satisfying 



Eq. (2.3 1 in Theorem 1. In this physical example, we do not need Eq. (2.1 1 as it 



should be, because it is a gapped system without any i-dependent orthogonalization. 
In order to prove Eq. ( 1.10 ), we write the Hamiltonian [9] for this system (see Fig.[T]) 



if2t , 



^ _ /lOJBoiir / COS a e ""sma 
'^^ ^ ~ 2 V e*"* sin a - cos a 



and its normalized eigenspinors [5] 



cos(q!/2) 
e**^" sin(a/2) 



'Pupi't) = ( „^^2^■^/^/0^ ) and <^down(t) 



e-*"* sin(a/2) 
— cos(q!/2) 



Their corresponding eigenvalues are /iajBohr/2 and — /iWBohr/2, respectively. From 
Eq. (1.3), the exact solution 



C0s(Ai/2) — [l(wBohr ^ n)/A\ sin(At/2) 
os(Ai/2) - [i(a;Bohr + f^)/A] sin(At/2) 



cos(a/2)e-*"*/2 
sin(a/2)e'"*/2 
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Proof for Eq. (1.10): Wc now evaluate 
(3.2) 



<Pdo 



(t) Hit) (^up(i)) = ^ 



■ sm a 



fku 



Bohre 



Since WBohr is a constant and f2 = 1 /To; 
given in Eq. (1.10) ■ 



we can readily obtain the proportionality 



If cJBohr is i-dependent, then one obtains 

•i^down 

(i) Lff(t)U„p(i) 



(3.3) 



Enp{t) — i?down(i) 



Equation (3.3 1 implies that Tex 3> 1 for QAA to be valid, and this inequality does 
not refer to Tin explicitly as explained in the introduction. It is not explicit because 
we can rewrite Eq. (|3.3|) to obtain 

h 



(3.4) 



,(t) Hit) ^„p(t))oc^/^oc 



« 1, 



Te.it) TUt)TeAt) 

which means that we are still not able to find the relationship between Tmit) and 
Tex(i). Therefore, Eq. (2.3) needs to be invoked independent of Eq. (2.2). 

Example 2: Spin-^ tv^ro-level oscillating system. In the above example, 
flit) = Qt, whereas a and WBohr are ^-independent constants. On the other hand, 
the oscillating Schwinger Hamiltonian [Hi [27] 

^Bohr(<) f cosa(i) e-'"(*) sina(t)^ 
s*"(*) sina(i) 



Hsit) 



cos ait) 

and n, a and WBohr are all time-dependent variables. We again assume the two 
energy levels do not cross and there is no i-dependent orthogonalization. Here, 
Fig. [2] schematically captures the physical system represented by Hsit). In this 
case, we have ait) = Q!(i)i, a(t)2,- ■ • ; ct = di,- • • ; 17(i) = f7(t)i,- ■ ■ ; fl — fli,- ■ ■ ; 
'^Bohr(^) — i^Bohr(i)i,- ' ' ; '^Bohr(i) — '^Bohr(Oir ' ' ■ This mcaus that ouc will obtain 
many a(t), a, ^it), 51, WBohr(i), i^BohrC^), and so on due to t-dependent fluctuations 
in the rotating magnetic field (see Fig. [2|. Now, to guarantee QAA for Hgit), we 



just need to invoke Eq. (2.3) such that WBohr(Oa: ^ \'^x\, |'^Bohr(i)ai 
and a;Bohr(i)« > ^u, \i^Bohrit)v\ > \^v\, ■■■ where {x,y,u,v} = N*. 
inequalities, we prove that the two conditions 



> \a. 



Using these 



(3.5) 
(3.6) 



J7 sin a — ia 






(l cos a — C^Bohr ~ Jt^'"S(^ ^^"^ ^ *"^) 







^ 1, and 





d n' 


/' 

'0 





dt! < 1, 



derived by Comparat [TT] are actually the extended versions of Eq. (2.3) 
Proof: For Holt), Eq. ( [3!5| reduces to 

fl sin a 



(3.7) 



cos ( 



~ l^Bohr 

ilt and 



«1, 



CI. Equation (3.7) requires 



because d = = (d/dt)arg(fl'), ^(t) 
WBohr ^ ^ for QAA to be valid, in accordance with Eq. (2.3). We now leave 
Eq. (3.5) intact and for QAA to be valid in Hsit), one requires wboIu- ^ |d|, WBohr ^ 
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Figure 2. A constant magnetic field, B (magnitude Bo) rotates 
with some fluctuations at a velocity, D,{t) and with a t-dependent 
opening angle, a{t). See Example 2 for details. 



n and |wBohr| ^ \^\ in which, the first two inequalities satisfy Tj^ <^ Tpx, while 
the last one is equivalent to iTlnl < jTexl- If (d/dt)arg(ri') < in Eq. (3^), then 
one further requires |a;Bohr| ^ |(d^/di^)arg(ri')| or jTinl <^ |Tex|, which guarantees 



WBohr ^ I(d/di)arg(r2')| as required from Eq. (3.5). Having proven that, it is now 



straightforward to generalize Eq. (3.6) to obtain 



(3.8) 



T 



< 1. 



Therefore, Eqs. (3.5) and (3.6 1 are indeed the extended versions of Eq. (2.3) 



Example 3: Two-level system v^rith a single energy level crossing (logical 
exposition). For simplicity, we confine ourselves to a two-level system with one 
energy level crossing occurred somewhere between to and ti, say tz- Again, we 
assume there is no i-dependent orthogonalization. In this case, the energy level 
crossing occurs as a result of some i-dependent Yukalov (non-linear)-type interac- 
tion, V{t) [TT. Hence, one can write the solved Hamiltonian for such a system 
as 



(3.9) 



H{t)^^{t) ^ [Hoit) + V{t))iprn{t) ^ {ha{t) + Va{t))ip,n{t) = Eait)ipm{t), 



where Ho{t) is the basic Hamiltonian excluding all interactions, whereas V{t) con- 
tains all the interaction terms, and we do not require {fn{t)\Vmit)) = or 1. 
Here, Ea{t), ha{t), Va{t), hb{t), Vb{t) and Eb{t) are all eigenvalues. Degeneracy 
due to energy level crossing imphes Ea{tz) — Eb[tz) at a certain point of time 
{tz), and for other times, they are not degenerate. If they are always degener- 
ate, then ha{t) = hb{t), Va{t) = Vb{t) and Ea{t) = Eb{t), which physically mean 
ipn (t) = (t) ■ this means that there are two electrons occupying the same energy 
level for all time t. As a consequence, degeneracy due to energy level crossing re- 
quires ha{tz) 7^ hb{tz), Va{tz) 7^ Vb{tz) and Ea{tz) = Eb{tz). Therefore, for degener- 
ate systems, one needs to employ Eq. (2.3) such that 1/Tex(tz) ^ {ha{tz) — hb{tz)) /f^ 
and 1/Tcx(tz) <C {va{tz)-Vb{tz))/n^\iem {va{tz)-Vb{tz))/h and [ha{tz)-hb{tz)) /h 
equal l/Ti-a^{tz)- This is the reason why QAA can still be apphed for gapless (due 
to energy level crossings) systems without strictly requiring Tex — ^ oo. 
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Example 4: Systems with and without the time-dependent orthogonal- 
ization: applied to chemical reactions. In the previous examples, we have 
shown how and why QAA is valid for gapped and gapless systems provided that 
they satisfy Eqs. (2.2) and (2.3). These conditions are also shown to be necessary 
and sufficient. In the following example however, we will first provide the logical 



expositions of how Eq. (2.1) due to i-dependent orthogonalization can render the 



QAA valid. All we need to do here is to invoice Eq. (2.1 1, which now implies 



(3.10) 



fdovjn{tl 



Hit) ^,p(to) 



Eup{to) — -E'down(il) 



1. 



Equation (3.10) tells us that the system has started out with spin-up, (fupito) att — 
to, and after t-dependent orthogonalization, it ends up with spin-down, (/7down(ii) at 
t = ti where ti > to. For QAA to be valid, we require large external disturbances 
such that T^^{t) w 1 or T^^{t) ~ T(^{t), which is now necessary for the vaHdity 
of QAA. Without this large disturbances (or rapidly oscillating Hamiltonian) , the 
spin- up electron (at t — to) cannot occupy the orthogonalized spin-down state at 
t — ti. In this case, c{ti) ^ 1 (invalid QAA) as proven in Eq. (2.20) and explained 
in the paragraph after Eq. (2.20). 



Claim: However, one should note here that Eq. (2.1 ) is also applicable for systems 
without the t-dependent orthogonalization. 
Proof: Recall Eq. ( [216] ) with St.to = 1 

Cm(il) = "Cm(ii) {iprniti)\iprn{tl)) - 



(3.11) 



5Z 



Xto] 



tpniitl) H (Prnito)) - E„i{to) 



„i[e„(to)-e„(ti)] 



E,n{to) — Em{tl) 



Here, St-^to = 1 implies that there is no i-dependent orthogonalization. First we 
neglect the second term by requiring {ipm(ti)\H\ipfn{tQ)) — i?,„(to) is extremely 
small and/or Em{to) — E„i{ti) is extremely large, hence 



(3.12) 

Its solution 
(3.13) 



Cm{tl) = -C„i(ti) {ipm{tl)\iPm{tl)) 



Cm(ii) = c„(to)e''^'"^*'\ where 



7m(ii) = i 



fm{t'] 



dt 



-ip,n{t'))dt' 



Since the electron occupied the m**^ eigenstate at t = to, thus Cm{to) — 1, and when 
t = ti, Cm(ti) = 1, and so on. Therefore, we can now recall Eq. (|2.5[) and write 



(3.14) 



^rn{to) ^ C^ita)^r,^ito)e''^"^^'°\ 



and using Eq. (3.13) 

(3.15) = c„(ii)^,„(ti)e'^"(*i) = c™(to)^™(ii)e^^™(*^)e'^"(*i). 



Since Stj^tg = Eq. (3.15) reads 

(3.16) vl/^(ti) = cMvMe''""^''^e''"'^''\ 



14 



ANDREW DAS ARULSAMY 



as it should be, picking up only an additional phase factor (e*'''™(*i)). If the second 
term on the right-hand side of Eq. (3.111 is not extremely small, then we strictly 
require 



(3.17) 



LP, 



Em{to) - 



Em{ii. 



Emita) 



oc 



1 



in accordance with Eq. (3.101 



We anticipate that the evolution of the eigenfunction in the presence of, or in the 
absence of t-dependent orthogonalization is relevant to understand chemical reac- 
tions beyond the standard procedures of calculating activation energies, including 
the molecule ion dissociation [TH] and the water molecule splitting [55]. For 
example, chemical reactions (associations or dissociations) strictly require large ex- 
ternal disturbances that give rise to significant changes to the eigenfunctions to 
form new compounds or to dissociate to its constituent chemical components |29j . 
In other words, without the significant changes to the eigenfunctions due to large 
external disturbances, which define the existence of activation energies, the chem- 
ical reactions cannot occur. For instance, bringing an atomic hydrogen and 11+ 
together chemically to produce requires a new wave function, if{t, r)ncw for 
the ground state of molecule. We take the wave function for the atomic H 



as (po{t,r). Thus, the condition given in Eq. (2.1) considers the evolution of the 
wave function, from ipQ(t,r) to (p(t,r)„ow However, the generator (Q) defined by 
ip{t, r)ncw = G'i^oit, r) is not known, and presently Q stands for guess as we often do 
exactly that to obtain the new wave function by means of the linear combination of 
atomic orbitals and/or by requiring any arbitrary wave function to give solutions 
that are convergent. Anyway, we can now invoke Eq. (2.1), which is also suitable 
to check the applicability of the QAT for the chemical reaction, H -I- 11+ — >■ H J . In 



this case, Eq. (2.1) can be written as 



iVchcmical 
reaction 



r 

E,n(il) — i?m(Hj) 



1, 



(3.18) 

where ^^^(H) and _Em(Hj) are the ground state energies of atomic H and mole- 
cule, respectively. The magnitude of approximately one stated in Eq. ( |3.18[ ) can be 
achieved during chemical reactions because large external disturbances (by means 
of thermal and/or potential energies) are required to initiate any chemical reaction. 
Such a disturbance will first give rise to excitations of the ground state electrons 
and therefore, QAA is also applicable for such a dynamical process, provided that 
we know the Hamiltonian (^f(i)rcaction') describing this chemical process explicitly. 



On the contrary, Eq. (2.2) is not applicable for this process, and as we have stated 



earlier in the introduction, this does not imply that the QAT has been violated just 
because ip{t,r)acw diff'ers significantly (beyond the phase factor) from (pQ{t,r). 

4. Further Analysis 
In the discussion and proofs given above, we should not be falsely led into think- 



ing that the condition given in Eq. (2.1) imposes orthogonality. On the contrary. 



Eq. (2.1 ) allows the final eigenfunction to remain the same, or acquire some phase 



factors or be mathematically different from the initial eigenfunction, or even be 
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different such that the final cigcnfunction is orthogonal to the initial one. More- 
over, the proof for Eq. (2.1 1 also correctly exposes the intrinsic reason why we 
need to guess the eigenfunction, for instance, we have used the linear combination 
of the atomic eigenfunction in the form of Eq. (2.5). In other words, we need to 
understand what is the mechanism involved, which gives us the license to rewrite 
(po{t,r) as ip{t,r)ncw Well, of course ip{t, r)ncw is one of the acceptable solution to 
the Hamiltonian that gives the lowest ground state energy. But we still lack the 
fundamental knowledge that enabled us to construct ip{t,r)^cw by hand. 

For example, we reconsider a molecular ion such that the single electron 
molecular orbital, ip{t, r)ncw can be constructed, either from the molecular-orbital 
(MO), or valence-bond (VB) theory [29 . In MO theory one uses the linear com- 
bination of atomic orbitals (LCAO) approach to construct the molecular orbital, 
namely, (y9(i,r)„ew is constructed from the atomic orbitals of the type, ipo{t,r). 
In this case, one immediately starts the construction in the molecular framework, 
meaning, in the presence of both H and 11+ . In VB theory, one still uses the atomic 
orbitals, but they are combined non- linearly [35]. Therefore, in both cases, the 
constructions of ip{t^ r)nQw are carried out in the molecular framework [29| . even 
though they are generated from the atomic orbitals in different ways. 

Now, the chemical reaction between H and H+, physico-chemically means that 
we are bringing two isolated atomic H and H+ ion closer from infinity, or they may 
come closer on their own due to some attraction between them via the Coulomb 
interaction at finite distances. Either way, once they are close enough, and with 
sufficient external heat (Q), but not higher (than the H— 11+ bond energy), the 
single electron in H can be excited (or polarized) to the highest possible atomic 
energy level and binds to the H+ ion, to form H^. When one reads the above 
chemical reaction carefully, two intrinsic phenomena pop out that can be written 
as two statements. 

Statement 1: There is a change in the mathematical property of the wave 
function (beyond the phase factor) when ipQ(t,r) for H transforms to ip{t,r)ncw for 

Remark 4 ■ Let us call this the wave function transformation, and this transfor- 
mation is one form of electronic phase transition arising due to the transition from 
the atomic (or ionic) to molecular system. The formal proof for the existence of 
such a quantum phase transition in many-body systems is given elsewhere |30j . 

Statement 2: In both MO and VB theories, ip{t, r)ncw have been constructed 
based on the assumption that the chemical reaction has already taken place, since 
the construction of wave functions were carried out within the molecular framework. 

Remark 5: This means that, we are lack of the information needed to find the 
generator (Q) pointed out earlier, and as anticipated, Q still remains as a guess, 
while to those who are not comfortable with this, we may call it "an educated 
guess" . 

In this work, we have actually proven these two statements by proving the ex- 
istence of the condition stated in Eq. (2.1). For example, we can observe that 
Eq. (2.1) correctly captures the effect of the chemical reaction by requiring an 
electronic phase transition (due to wave function transformation), in which, the 
existence of this transition justifies why Remark 2 is necessary. Simply put, even 
though time is a continuous parameter, we cannot construct a continuous t- or 
r(i)-dependent eigenfunction for an electron that starts as an atomic orbital of an 
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isolated atomic H that goes through an electronic phase transition to form HJ. 
Here, r{t) denotes the t-dependent electron coordinate, which means, we do not 
require ip{t, r)ncw to be both r and t dependents explicitly. 

As a consequence, just because ip(t, r)ncw cannot be written as a continuous 
function of r{t) (starting from isolated H and H+, imtil the formation of H^), one 
should not assume that the above-stated wave function transformation (or phase 
transition) never took place, or should be ignored. The very moment we combine 
the atomic orbitals (linearly or otherwise) within the molecular framework proves 



the existence of this transformation, which we discovered from Eq. (2.1 ). Moreover, 
Eq. ( |2.1[ ) also leads us to understand why the activation energy for the above chemi- 
cal reaction exists— because it is needed to initiate the wave function transformation 
(non-observable) or the electronic phase transition (observable). This is consistent 
with the finite-temperature quantum phase transitions proven in Ref . |30j . 

Finally, before we conclude, let us recall the Born-Oppenheimer (BO) arguments 



given by Teufel [5|. Here, Eq. (2.1 1 does not apply for an obvious reason— a given 
wave function in a particular quantum system, that does not go through a quantum 
phase transition or f-dependent orthogonalization, does not transform beyond the 



phase factor with time. Therefore, we just need to invoke Eq. (2.2) and Eq. (2.3) 
for such systems, which can be done within these well-known theories, MO, VB 
and the Density Functional Theory (DFT). Briefly, BO approximation is valid 
for as long as the electrons average kinetic energy equals, or larger than that of 
the ions, ■^mv'^i > ^Mvf^^, which in turn implies Vion ^ Vc\ and Tox ^ ?in if 
M ^ m. Obviously, the above assumption does not violate Eq. ( |2.3[ ), hence QAA 
is applicable. 

It is also true that repulsive interaction (due to electron-electron interaction) can 
be invoked to explain the existence of finite energy gap [51 [3T], however, energy levels 
do cross due to "lucky coincidences" in many-body systems, as shown in Example 
3. Even in the presence of such coincidences, we can still check the applicability of 



QAA via Eq. (2.3). Now, in the worst case scenario, we will have (refer to Example 
3) 

(4.1) hait) = hb{t), Va{t) = Vb{t) and Ea{t) = Eb{t), 



then QAA seems to be not applicable. But if we look closely, Eq. (4.1 1 physically 



imply (/?„ — iprn- Hence, one can indeed invoke Eq. (2.3) through the BO approxi- 
mation and justify the applicability of QAA because the transition between Lpn and 
ifm does not violate QAA. Of course, mathematically one can construct the eigen- 



functions such that 7^ </j„i and still satisfy Eq. (4.1). But this is like saying, we 
are attempting to invalidate QAA, solely for the purpose of violating QAA without 
any physical justification. 



5. Conclusions 

In conclusion, we have developed a new version of the quantum adiabatic the- 
orem that takes the time-dependent orthogonalization into account, and even in 
the absence of this orthogonalization, we found that these (with and without the 
^-dependent orthogonalization) effects are very relevant to study the dynamics of 
chemical reactions. We also generalized and incorporated the relationship between 
the internal and external timescales as an additional, but separate condition into 
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the quantum adiabatic theorem. The time-dependent orthogonahzation and chem- 
ical reactions are proven to be captured by Eq. (2.1) self-consistently. In addition, 
Eq. (2.1) also led us to discover the activation energy in a given chemical reac- 
tion is needed for the wave function transformation, which can be observed via the 
electronic phase transition. While the validity of other conditions are explicitly 
exposed to be sufficient in the spin-1/2 two-level system. Therefore, the new the- 
orem is shown to be valid for both gapped and gapless systems with appropriate 
examples. 

For gapped two-level system, Eqs. (2.2) and (2.3) are found to be sufhcient and 
necessary, in accordance with Comparat (llj . Whereas, for gapless systems with 
degeneracies (due to energy level crossings) and quantized energy levels, Eqs. (2.2) 
and (2.3) are again sufficient and necessary, in agreement with Avron and El- 
gart [12 . However, for gapless systems, we need not keep track of the eigenfunc- 
tions to show they are sufficient, for as long as Eq. (2.3) is valid. Interestingly 
however, we do not always require Tox — >■ oo in the presence of non-linear Yukalov- 
type [17] interactions. The importance of Eq. (2.3) cannot be overstated as the 
recent counter-examples [HI |32] and counter-proofs |lT| [TH [32] can be traced back 
to the violation of Eq. [55] , 
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